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D ■ Abstract. In line with some previous works, we study in this paper the meson spectrum in 

' the framework of a second order quark-antiquark Bethe-Salpeter formaUsm which includes con- 

finement. An analytic one loop running coupling constant as{Q), as proposed by Shirkov and 
Sovlovtsov, is used in the calculations. As for the quark masses, the case of a purely phenomeno- 
logical running mass for the light quarks in terms of the c. m. momentum is further investigated. 
, Alternatively a more fundamental expression mp{Q) is introduced for light and strange quarks, 

^ ' combining renormalization group and analyticity requirements with an approximate solution of the 

\ Dyson-Schwinger equation. The use of such running coupling constant and masses turns out to be 

l/^ ■ essential for a correct reproduction of the the light pseudoscalar mesons. 

cn ■ 

(N : 

5^ : 1. INTRODUCTION 

O ; 

^ I In a series of papers [1, 2] we have applied a second order Bethe-Salpeter formal- 

ism \ previously established [3], to the evaluation of the quark-antiquark spectrum, in 
Qh! the context of QCD. Taking advantage of a Feynman-Schwinger representation for the 

r-| ! quark propagator in an external field, the kernels of the Bethe-Salpeter and the Dyson- 

■ Schwinger equations were obtained, starting from an appropriate ansatz on the Wilson 

•T-j . loop correlator. Such an ansatz consisted in adding an area term to the lowest pertur- 

\ bative expression of \nW . By a 3D reduction of the original 4D BS equation a mass 

' operator was obtained and applied to the determination of the qq bound states [4] . 

In that way, using a fixed strong coupling constant and appropriate values for 
the other variables, the entire spectrum was reasonably well reproduced with, however, 
the relevant exception of the light pseudoscalar mesons (tt,^, r\s). Agreement even for 
the latter states could be obtained using an analytic running coupling constant (Xs,{Q) 
proposed by Shirkov and Sovlovtsov, which is modified in the infrared region with 
respect to the ordinary purely perturbative expression [5, 6]. In conjunction it was also 
necessary to use a phenomenological running constituent mass for the light quarks u and 
d, written as a polynomial in the center of the mass quark momentum k [2]. 
In this paper we reconsider and improve the above procedure from two aspects: 
a) we evaluate the hyperfine — ^Sq separation for the light pseudoscalar mesons to 
the second rather than to the first order perturbation theory. 



^ second order in the sense of the differential equations 



b) we use running constituent masses for u,d and s quarks, obtained by an approxi- 
mate solution of the appropriate DS equations and analytic running current masses. 

As a consequence of a) a significant improvement is obtained in agreement with the 
data already with a phenomenological running mass for the light quarks. As for case b), 
preliminary calculations seem to provide results numerically similar to the above ones, 
but more satisfactory from the conceptual point of view. 

The plan of the remaining part of the paper is as follows. In Sect. 2 we briefly recall the 
second order BS formalism to establish notations. In Sect. 3 we discuss the DS equation 
and the 3D reduction of the BS equation. In Sect. 4 we consider the infrared behavior of 
the running coupling constant and obtain the corresponding running masses. In sect. 5 
and 6 we report our results and draw some conclusions. 



2. SECOND ORDER BETHE-SALPETER FORMALISM 

In the QCD framework a second order four point quark-antiquark function and a full 
quark propagator can be defined as 

H^'^\xuX2;yi,y2) = -^Trcoior{^i{xi,yi;A)A2{y2,X2;A)) (1) 

and 

H^^\x-y) = -^Trcoior(A(x,3; : A)) , (2) 

where 

(/[A]) = |mM^[A]e'^G[^]/[A], (3) 

Mf[A] = Detnj^^[l + gY^An{iYj djv — mj)~^ and A{x,y;A) is the second order quark 
propagator in an external gauge field. 

The quantity A is defined by the second order differential equation 

(D^D^+m2- iga^^F^v)A(^,y;A) = -5\x-y), (4) 

(<7^^ = § [7^, 7^] and Dfj^ = df^+ igA^) and it is related to the corresponding first order 
propagator by S{x,y;A) = {iy"^ Dv + m)A{x,y;A) . 

The advantage of considering second order quantities is that the spin terms are more 
clearly separated and it is possible to write for A a generalized Feynman-Schwinger 
representation, /. e. to solve eq. (4) in terms of a quark path integral [3, 1]. Using the 
latter in (1) or (2) a similar representation can be obtained for ^nd _ 

The interesting aspect of this representation is that the gauge field appears only 
through a Wilson line correlator W. In the limit X2 ^ xi, y2 ^ yi or 3; — > x the Wilson 
lines close in a single Wilson loop T and if T stays on a plane, ilnW can be written in a 
first approximation as the sum of its lowest perturbative expression and an area term 

ilnW^ ^g^jdz^ ^dz^'D^yiz-z) + (5) 



./..»/W-.-)|z-.'|/'.A{l-lA|i + (,-X)§]^}' 



The area term here is written as the algebraic sum of successive equal time strips and 
dz±_ — dz— {dz ■ r)r/r^ denotes the transversal component of dz. The basic assumption 
now is that in the center of mass frame (5) remains a good approximation even in the 
general case, /. e. for non flat curves and when X2y^xi,y2^yi ory ^x. 

Then, by appropriate manipulations on the resulting expressions, an inhomogeneous 
Bethe-Salpeter equation for the 4-point function H^'^\xi,X2;yi,y2) and a Dyson- 
Schwinger equation for H^'^'> (x — y) can be derived in a kind of generalized ladder and 
rainbow approximation. This should appear plausible, even from the point of view of 
graph resummation, for the analogy between the perturbative and the confinement terms 
in (5). We may refer to such terms as a gluon exchange and a string connection. 

In momentum representation, the corresponding homogeneous BS -equation becomes 

^ .,s f d^u t f, 1„ k^-u 1 k + u\ 

Af^ O'M") (-5''+*) ■ 

where we have set = 1; a, I? = 0, /iv; the center of mass frame has to be understood, 
P — (mfi, 0); ^p{k) denotes an appropriate second order wave function ^. 

Similarly, in terms of the irreducible self-energy, defined by H^'^'> {k) — j^il^i + 
1^2 1^2 i^{k) H^'^^ {k) , the DS-equation can be written 

m=j^L{^-,:'-±L,'-±Lymi)-^. (7) 

The kernels in (6) and (7) are the same in the two equations, consistently with the 
requirement of chiral symmetry limit [7], and are given by 

IodQ;p,p') = l6K^(Xsp"p'PDap{Q) + 

+4(7 Jd'Ce-''^<\C\e{po)e{p'o) l'dX{plp'i-[Xp'oPT+{l-^)poP^?y^ 
I^v;o{Q;p,p') = 4ni^as{d^Qy - d^Q^)p'^b^p {Q) - 

Po 



ICIVpo-Pt 



In terms of the second order field (j) (x) = (iY^Dj^ +m) ^ \i/{x) this wave function is defined by 



^ \^\^Jp'i-v'^ 



(8) 



where in the second and in the third equation Co = has to be understood. Notice that, 
due to the privileged role given to the c. m. frame, the terms proportional to <7 in (8) are 
not formally covariant. 



3. DS EQUATION AND MASS OPERATOR 

Concerning eq. (7), let us observe that the unity matrix, cr^^ and form a subalgebra 
of the Dirac algebra. Consequently r(^) can be assumed to depend only on this set of 
matrices and, since it must be a three dimensional scalar, only on terms like kjo'^K In 
fact, it can be checked that r(^) can be consistently assumed to be completely spin 
independent and eq. (7) can be written in the form 



dH R{k,l) 



(9) 



with 



R{k,l) 



, 4 



yi'-ipiy 1 3 

{k-lf 4 



+ 



+ 



'3rg-i(k-l)T 



•v(^o+/orwi 



(kj. + L)2 



(10) 



and denoting as above the transversal part of k and 1. 
Notice that, once (10) is solved, the pole or constituent mass mp, to be used in bound 
states problems, is given by the equation 



mp-m2 + r(mp) = 0. 
We can try to solve eq. (9) iteratively and we have at the first step 

dH R{k,l) 



m=ij 



(11) 



(12) 



In a preliminary calculation we omit altogether the perturbative contribution to R{k, I) 
(notice the overplacing of curves b and c in Fig. 1) and neglect the term in (k^ + 1^)^ 
in the string part. Strictly, the second approximation is justified only for S bound states 
(classically k±r is the angular momentum of the bound state) but it is necessary in order 
to make the integral analytically calculable. 



Then introducing a cut off jU, we obtain 



r{k) = -[klA{m, \k\)-B^-B{m, |k|) , (13) 

7t 



where5„ =21n^-l. 



A(m,|k|)= ^ 



+m 



2 



, Vk^ + + Ikl 

1 H In — ^ 

2|k|Vk2+m2 Vk2 + m2-|k| 



(14) 



and B{m, |k|) is a more complicated expression that we do not report explicitly here for 
lack of space. The resulting pole mass is 



m2 + f [5^+5(m, 


k|) 




k|)] 


1 + f A(m, 


k|) 



mm,\k\) = E^JL^^1LJ1^^^±J^ ^ (15) 



The above expression depends on the current mass m and on the quark c. m. momen- 
tum |k|, (see Fig. la). Notice that such dependence on |k| is clearly an artifact of the 
schematic way we have introduced confinement in eq. (5) and that the curve is rather 
flat in the region of interest. Correspondingly it seems reasonable to chose as true mass 
nipim) the value of mp{m, |k|) at its stationary point in |k|. 

Then, in a neighbor of its singularity — mj,, the full propagator can be written as 
^(2) (/j) — -j^^^, where the residuum Z differs from 1 only for terms proportional to Ug 

or o. Consistently in (6) we can simply take Z = 1 and are left with the free propagator 
with a constituent mass. If, in addition, we replace lab with its so called instantaneous 
approximation 4^^'(k, u), we can explicitly perform the integration in uq and arrive at a 
three dimensional reduced equation. 

Such a reduced equation takes the form of the eigenvalue equation for a squared mass 

operator [3], = Mq + f/ , with Mq^w\+W2, w\;i = W 2 + and 



(k|t/|k)-^y^^/„, (k,k)y^-^CTiCT2 (16) 

(for an explicit expression we refer to [2, 1]). The quadratic form of the above equation 
obviously derives from the second order formalism we have used. 

Alternatively, in more usual terms, one can look for the eigenvalue of the mass 
operator or center of mass Hamiltonian Hqm = M = Mq + V with V defined by MqV + 
VMq + = U. Neglecting the term the linear form potential V can be obtained 

from U by the kinematic replacement \ (>*'i+'^2)(vvi+w2) _^ — i j^iq resulting 

expression is particularly useful for a comparison with models based on potential. In 
particular, in the static limit V reduces to the Cornell potential 

Vstat = -^- + ar; (17) 
3 r 



in the semirelativistic limit (up to ^ terms after an appropriate Foldy-Wouthuysen 
transformation) it equals the potential discussed in ref. [8], if full relativistic kinematics 
is kept, but the spin dependent terms are neglected, it becomes identical to the potential 
of the relativistic flux tube model [3]. 



4. RUNNING COUPLING CONSTANT AND MASSES 

As we said, diagonalizing or Hqm with fixed coupling constant and quark masses, 
a general good fit of the data was obtained. Actually a serious problem was represented 
by the masses of the light pseudo scalar mesons that turned out too large. The results 
obtained in ref. [2] suggest, however, that the situation can be greatly improved using an 
appropriate running coupling constant and running quark masses. 
At one loop, the running coupling constant is usually written 

with j8o = 1 1 — \Nf and A'f the number of 'active' quarks. However, the singularity occur- 
ring in such expression is an artifact of perturbation theory and it contradicts general an- 
alyticity properties, therefore the expression must be somewhat modified in the infrared 
region [9]. Notice that this is particularly important for the quark- antiquark bound state 
problem, where the variable is usually identified with the squared momentum trans- 
fer q2 = (k - k')^, which ranges typically from (0.1 GeV)^ to (IGeV)^ for different 
quark masses and states. 

The most naive modification of eq. (18) would consist in freezing 0Cs(2^) to a certain 
maximum value as decreases and in treating this value as a phenomenological 
parameter (truncation prescription). However, various more sensible proposals have 
been made on different bases [5, 6]. 

In particular Shirkov and Solovtsov [6] suggest to replace (18) with 

. X 4;r / 1 a2 \ 

This remains regular for = A^ and has a finite A independent limit (x^{Qi) = 4-n/po 
for 0. Eq. (19) is obtained assuming a dispersion relation for OCs{Q) with a cut for 
— oo < <0 and applying (18) to the evaluation of the spectral function. 
The running mass expression corresponding to (19) can be written in the form 

where in the MS scheme 7o = 8. Eq. (20) is obtained integrating the one loop renormal- 
ization group equation 

Q dm{Q) _ ^^a,{Q) 
m{Q) dQ ~ 4;r ' ^ ^ 



where (19) has been used and m denotes an integration constant. Notice that m{Q) is 
singular for 2 — > 0, contrary to Cts(2)- 

Finally, if we replace the running mass (20) in (15) after maximizing we obtain a 
running constituent mass mp{Q) of the type reported in Fig. 1 d that can be used together 
with the running coupling constant (19) in the expression of the operator (see Sec. 
3).^ 



5. CALCULATIONS AND RESULTS 

The calculations we report in this paper follow a similar line to those of Ref. [1]. The 
general strategy for solving the eigenvalue equation for and the numerical treatment 
are basically the same. 

We neglect spin-orbit terms, but include the hyperfine terms in U (see eq. (16)). We 
solve first the eigenvalue equation for Mstat — Mq-\- Vstat (see Eq. (17)) by the Rayleigh- 
Ritz method with an harmonic oscillator basis and then treat — M^^^^ as a perturbation 
(up to the first order this is obviously equivalent to taking m| = (M^)). 

In the above general framework, in Fig. 2 we graphically report and compare with the 
data [10] three different type of results, corresponding to different choices for the strong 
coupling constant ttg, the string tension o and the constituent masses. 

Diamonds correspond to results already reported in [1]. A running coupling constant 
Cfs(2) was assumed equal to the one loop perturbative expression (18) frozen at the 
maximum value tts = 0.35, with Nf — 4 and A = 200 MeV. Q was identified with |k — k' | 
and a was set equal to 0.2 GeV^. Fixed masses niu — = 10 MeV, = 200 MeV, 
nic = 1.394 GeV, = 4.763 GeV were adopted. The results do not differ essentially 
from the fixed coupling constant case; the spectrum is reasonably well reproduced on 
the whole with the exception of the light pseudoscalar mesons n, rjs and K (the ris mass 
is derived from the masses of T] and r]' with the usual assumptions). 

Circlets correspond to results of the type reported in [2], but in which the hyper- 
fine separation for the IS and 2S states has been evaluated up to the second order 
of perturbation theory. In this case as running coupling constant we have taken the 
Shirkov-Solovtsov expression (19). We have set again Nf = 4 and A = 200 MeV, but 
(7 = 0.18 GeV^ and = 0.39 GeV, nic = 1.545 GeV, = 4.898 GeV. On the con- 
trary for the light quarks we have taken a phenomenological running mass in terms of the 
modulus of the cm. quark momentum, = rn^i = 0.17|k| — 0.025|kp + 0.15|k|^ GeV^ 
with Ik I in GeV. We can see that in this way even the light pseudoscalar mesons turn out 
correctly, with possibly some problems for the qb states and some other highly excited 
states for which coupling with other channels are probably important. 

Squares correspond to preliminary completely new calculations, made using the an- 
alytic running constant (19), running constituent masses mp{Q) (as described in Sec. 



^ At first sight it could seem strange that we should talk of a 2 dependence for a quantity like the 
constituent mass that should have a definite physical value. The point is that we are using mp{Q) in 
the context of certain approximations and it is the accuracy of such approximations that depend on the 
scale Q. 




4) for the light and strange quarks, fixed masses for the charm and the beauty quarks. 
Inside (Xs{Q) the quantity Q has been again identified with |k — k'|. On the contrary, for 
computational difficulties, inside mp{Q) we have taken 

2=^, (22) 

where Ye is the Euler constant ye = 0.5772 . . . and (r) is the radius of the unperturbed 
bound state [11]. We have chosen % = 3, A = 180 MeV, o = O.lSGeV^, §5^ = 0.48 
GeV in (15), both for the light and the strange quarks, and then m„ = = 25.0 MeV, 

= 87.3 MeV in (20) (in order to reproduce correctly the p and the ^ masses). Finally 
we have taken = 1.508 GeV and = 4.842 for c and b quarks. The results are 
not of a better quality than those obtained in the preceding calculation but obviously 
conceptually more satisfactory. 

As an example, in the table numerical values for the three types of calculations are 
reported in the order for the light-light channel. For the third case in the last column the 
pertinent values of running constituent light quark mass are also reported for the various 
states. 



6. CONCLUSIONS 

In conclusion we can confirm what already noticed in references [2] that our reduced 
second order formalism together with ansatz (5) can reproduce reasonably well the 
general structure of the entire meson spectrum, light-light, heavy-heavy and light-heavy 
sectors included. In order to obtain the masses of light pseudo scalar mesons K, rjs and 
K, however, a correct consideration of the infrared behavior of the running coupling 
constant and of some kinds of running constituent mass for the light quarks is essential. 
The analytic Shirkov-Solovtsov coupling constant seems to provide such a behavior. 

What is new in this paper is the inclusion of second order perturbative corrections to 
the hyperfine splitting in the case of phenomenological running masses considered in [2] 
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FIGURE 2. Quarkonium spectrum (lines represent experimental data). 



(circlets in Fig. 2) and the use of a running mass obtained combining renormalization 
group and analyticity requirements with an approximate solution of the quark Dyson- 
Schwinger equation (squares in Fig. 2). 



ACKNOWLEDGMENTS 

We are indebted to C. Simolo for various calculations on the running masses. 



REFERENCES 

1. M. Baldicchi, G.M. Prosperi, Phys. Rev. D 62 (2000) 114024; Fizika B 8 (1999) 2, 251; Phys. Lett. B 
436 (1998) 145. 

2. M. Baldicchi and G. M. Prosperi, Phys. Rev. D 66 (2002) 074008; Color confinement and hadrons 
Quantum Chromodynamics, Page. 183, H. Suganuma, et al. eds. World Scientific 2004, hep- 
ph/0310213. 

3. N. Brambilla, E. Montaldi, G.M. Prosperi, Phys. Rev. D 54 (1996) 3506; G.M. Prosperi, Problems of 
Quantum Theory of Fields, Pag. 381, B.M. Barbashov, G.V. Efimov, A.V. Efremov Eds. JINR Dubna 
1999, hep-ph/9906237. 

4. For a different approach to the relativistic bound state problem see : Yu.S. Kalashnikova, A.V. 
Nefediev, Yu.A. Simonov, Phys. Rev. D 64 (2001) 014037; and references therein. 

5. A.C. Mattingly, P.M. Stevenson, Phys. Rev D 49 (1994) 437; S. J. Brodsky, hep-ph/0412101; S. J. 
Brodsky, S. Menke, C. Merino, Phys. Rev D 67 (2003), 055008; P Boucaud et al, JHEP 04 (2000) 
006; A. Ringwald and F. Schrempp Phys. Lett. B 459 (1999) 24; Yu.L. Dokshitzer, A. Lucenti, G. 
Marchesini, G.P Salam, JHEP9805 (1998) 003; Yu.L. Dokshitzer, V.A. Khoze, S.I. Troyan, Phys. Rev 
D 53 (1996) 89; N.V. Krasnikov, A.A. Pivovarov, Phys. Atom. Nucl. 64 (2001) 1500; G. Grunberg, 
Phys. Rev D 46 (1992) 2228; A.l. Sanda Phys. Rev Lett. 42 (1979) 1658. 

6. D.V. Shirkov, I.L. Solovtsov, Phys. Rev Lett. 79 (1997) 1209; Theor Math. Phys. 120 (1999) 1220; 
D.V. Shirkov, hep-ph/0408272; A.I. Alekseev, B.A. Arbuzov, Mod. Phys. Lett. A 13 (1998) 1747; A.I. 
Alekseev, Few Body Syst. 32 (2003) 193; see also A. V. Nesterenko, Int. Joum. Mod. Phys. A 18 
(2003) 5475. 

7. M.B. Hecht, CD. Roberts, S.M. Schmidt, Phys. Rev C 63 (2001) 025213, and references therein. 

8. A. BarchieUi, E. Montaldi, G.M. Prosperi, Nucl. Phys. B 296 (1988) 625; Erratum-ibid. B 303 (1988) 
752; A. BarchieUi, N. Brambilla, G.M. Prosperi, IlNuovo Cimento 103 A (1990) 59; N. Brambilla, P 
ConsoH, G.M. Prosperi, Phys. Rev D 50 (1994) 5878. 

9. I.F. Ginzburg, D.V. Shirkov, Sov Phys. JEPT22 (1966) 234. 

10. The Review of Particle Physics, S. Eidehnan et al., Phys. Lett. B 592 (2004) 1. 

11. W. Lucha, F. Schoberl, D. Gromes, Phys. Rep. 200 (1991) 127. 



